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Renewable energy based microgrids have main challenges of maintaining its 
frequency-voltage characteristics and system becomes more complex when 
they are interconnected. These sources being intermittent in nature need to be 
supported by other resources like diesel/biogas such that at time of small 
variation in load or natural sources (wind/solar), power requirement is met 
through support provided by diesel/biogas-based system. Also, the controller 
should be fast enough to minimize the changes such that system reaches 
steady state. In this paper, renewable based rural microgrid consisting of 
wind, solar and biogas is modeled and interconnected through low voltage 
AC (LVAC) line. Also, one of the microgrid modeled is connected to the 
main grid as well as drawing power from the other microgrid. Control 
approach have been developed in such a way that whenever there is 
disturbance in the system due to increase/decrease in load or input to the 
renewable energy sources the biogas-based system of individual microgrid 
increases/decreases its generation to support the system requirement. No 
extra power is drawn either from the LVAC network or main grid as desired. 


modeling of system and its dynamic Study has been carried out in 
MATLAB/Simulink. 
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1. INTRODUCTION 

Renewable energy generation-based power plants in small and large scale are needed to be installed 
in India to meet the power requirement and decrease dependency on thermal [1]. In this paper dynamic 
stability study has been carried for cluster of villages forming rural Microgrid with local load and 
interconnection with other rural microgrid as well as main grid. Novelty in the paper is modeling based on 
realistic approach in formation of microgrid based on available renewable energy sources and its control. 
Such type of microgrid is required in future to minimize long distance transmission, reduce dependency on 
main grid and to harness renewable energy sources in small scale. 

Use of power electronics has made it easier for asynchronous interconnection of microgrid through 
dc link and back to back convertor and also drawing power from WES/PVS [2]. The control and optimization 
of a hybrid ac/dc network are still an open problem [3]-[5]. The coordinated control architecture using power 
electronics devices has made it easier to control the system frequency and voltage [6], [7]. Inverters form 
critical equipment in control architecture of controlling power from renewables [8]. Many studies have been 
carried out in active and reactive power control of hybrid power system [9]-[11]. Therefore, generation mix 
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of wind-solar-biogas could be developed to provide power to habitants and small-scale industries locally 
using optimal use of renewable energy resources. Various literature and works have been carried out 
considering wind turbine with different types of generators and its control. The main goal of any system is to 
maintain a constant voltage profile and frequency produced by the system and this can be expressed in terms 
of constant angular velocity of the turbine rotor. Therefore, the angular position of the rotor works as a 
control signal. PMSG with wind turbine can be designed to generate high voltage at small revolution and 
high current at faster speeds and also has many advantages [12], [13]. Wind turbine with Induction Generator 
requires reactive power control which is supplied by capacitor bank or static compensator [14]. PVS dc 
power output dependents on temperature and solar irradiation. In order to excerpt extreme power from a PV 
cell at a given temperature and solar irradiance operating conditions, a maximum power point tracking 
(MPPT) system is essential. The most two common maximum power point tracking algorithms are: perturb 
and observe (P&O) and incremental conductance (IC) [15], [16]. 

Microgrid is energy system which is locally controlled and uses various categories of renewable 
energy assets (solar, wind, biomass, hydro and ocean), generators (biogas, diesel and gasoline), energy 
storage systems (batteries, flywheel, fuel cell and thermal), loads (residential and industrial) and control 
equipment [17]. A microgrid model optimizing local renewable energy resources for on-grid area was 
investigated for energy cost saving and reduction in CO» emission [18]. The problem of interconnecting 
microgrids is major area of concern because being small system, frequency, voltage and phase angle may 
vary with nominal disturbance. Since the microgrid consist of various generating sources, due to different 
power sharing capacity the chances of microgrid failure are higher in case of fault or heavy change in load or 
in drop in renewable energy source [19], [20]. Main cause of collapse is very fast response of inverter-based 
system and slow response of prime mover based reciprocating engine. In such case controller is designed and 
tuned in such a way so that effect is nullified very fast. Fractional order proportional integral controller 
(FOPI) (PI?) based controller is one of the solutions. More flexibility was extended in traditional notion of 
PID controllers with the controller gains define the fractional differential and integrals as design variables 
[21]-[23]. The fractional-order system can be approximated by many methods [24]. Based on approximation, 
the FOMCON toolbox was used for creating the fractional-order system [25], [26]. 

The rest of the paper is organized as follows: Section-2 gives the detail of the system considered for 
investigation. Section 3 comprises of description regarding modeling of microgrids and its associated line. 
Section 4 discuss about the simulation of the system and results. Finally, section-5 is the conclusion of the 
paper under consideration. 


2. SYSTEM CONSIDERED FOR INVESTIGATION 

In the system modeling two clusters viz. MG-1 and MG-2 are considered for stability study. System 
details with FOPI tuned parameters, low voltage AC (LVAC) line data, MG-1 data and MG-2 data are given 
in Table 1, Table 2, Table 3 and Table 4. Both the clusters are rated at 1500 kW. It is assumed that the power 
is flowing from MG-1 to MG-2. MG-1 consists of WES rated at 1000 kW and biogas generation system 
(ВОТ) with synchronous generator rated at 500 kW. Data of BG1 and WES are given in Table 5 and Table 6. 
МО-1 has local load of 1000 kW and also transmitting 100 kW to MG-2 considering loss of 5 kW through 
LVAC network i.e total generation of MG-1 is 1105 kW. MG-2 consists of PVS rated at 650 kW, Biogas 
generation system (BG2) rated at 500 kW and grid rated at 350 kW. Data of BG2, PVS and grid are given in 
Table 7, Table 8 and Table 9. LVAC network consists of 3-phase overhead Ariel Bunched cable (ABC). The 
size of cable considered is 150mm? aluminum conductor rated at 415 Volts. The objective of the LVAC 
network is to supply constant power of 100 kW to MG-2. Supply from limited grid is considered constant 
and is supplying to priority loads only. In the limited grid the voltage as well as the angle at that bus is kept 
constant. 


Table 1. FOPI data 


Controller мса Mo? 
P(Ky) KK) a P(K,) KK) a 
FOPII 6.987 2.578 0.45 2.745 4.572 0.8 
FOPD -6.296 -8.578 05 -1.425 -7.887 0.71 
FOPI3 1.945 10.676 0.8 10.355 2.765 0.2 
FOPI4 -2.378 | -5647 0.7 9.567 -3097 0.9 


The MG-1 and MG-2 microgrid has been modeled considering total contract demand=1725 kW. 
Maximum differentiated demand at diversity factor of 0.58=1725*0.58=1000 kW (approx.). 
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Table 2. LVAC data 
value value value value 
R=0.562 pu X.0.268 pu K jac=1.758 Kyac=-1.137 
Кз«=1.132 K44,-0.252 К5ж«=-1.132 K6ac=-1.137 
Kyac=1.514 Kga-=-1.381 Koac=0.531 K joac=0.531 
Ky tac=- 1.733 K joac=- 1.514 
Table 3. Generation parameters of MG-1 
Sources Rated Capacity (kW) Generation (kW) Participation Factor 

Biogas 500 350 0.2333 

WES 1000 755 0.5020 

Total 1500 1105 


Drs! = P1 0.01333 pu kW/Hz. 


Of Prf 


Krsi-1/Dgg1775 Hz/pu kW. 


Trsi-18 sec. 


д 
Dyn ecu E 


OV Орки Pg — 


0.6777 pu КУАВ/ри КУ. 


Kvsi-1/Dy5,71.5 pu КУ/ри КУАК. 
Туѕ-0.004 sec. pf=0.866, PL 1=0.6667pu 


Table 4. Generation parameters of MG-2 


Sources Rated Capacity (kW) Generation (kW) Participation Factor 
Grid 350 300 0.2 
Biogas 500 300 0.2 
PV 650 300 0.2 
Total 1500 900 +100 (Line) 
Desas = 77120,0133 pu kW/Hz. 
m 
Kro-1/Drsz =75 Hz/pu kW. 
Trs2-6 sec. 
дог Qin 2 Рр 


М: 


ôV Орки PgsV 


—0.6667 pu kVAR/pu kV. 


Kyo-1/Dys2=1.5 pu kV/pu КУАК. 
Ty2-0.004 sec. pf=0.866, Pj ?20.6667pu 


Table 5. Biogas (BG1) data 


Table 6. WES data 


value value value value 
Рвс:=0.2333 рч kW Овс!=0.1346 ри kVAR Pwes=0.503 pu Qwes=0.290 pu 
Ха=1 pu Ха’=0.15 pu Химу=0.5 pu Т.=0.005 sec 
Tpoi=0.01 sec Твс2=0.02 sec Kwes1=2.265 Kwes2=2.265 
Твсѕ=0.014 5ес Твсв=0.04 ѕес Kwes3=0.433 Kwes4=0.5 
Tpo7=0.036 sec Kap=200 Kwess=-0.5 Kweso=-0.5 
Тьв=2.0 sec Кьв=0.5 Kwes7=1.953 Kwess=- 1 .735 
Tap=0.05 sec Кев=1 
Тев=1.0 sec Kipo1=0.15 
К›вс1=0.831 К;звс1=6.515 
Kypoi1=-6.844 Тао‘ =5.0 5ес 


Тсвс!=0.75 sec 


Table 7. Biogas (BG2) data 


Table 8. PVS data 


value value value value 
Pgoi1=0.2 pu kW Qs c170.1 154 pu kVAR Ppyi=0.2 pu Qpyi=0.1 154 pu 
Ха=1 pu Ха’=0.15 pu Хтру=0.8 ри Tpyi=0.005 sec 
Tpoi=0.01 sec Твс=0.02 ѕес Kpy;=1.467 Kpy2=1.467 
Твс5=0.014 sec Твсв=0.04 ѕес Kpy3=0. 169 Kpy4=0.2 
Tpo7=0.036 sec Kap=200 Kpys=-0.2 Kpyo=-0.2 
Trp=2.0 sec Krp=0.5 Kpy7=1 238 Kpys=- 1.03 
Tap=0.05 sec Кев=1 
Trp=1.0 sec Kipoi=0. 15 
К›вс1=0.831 К;звс1=6.515 
Kypoi=-6.844 Тао‘ =5.0 5ес 


Тсвс!=0.75 sec 


Table 9. Grid data 


value value 
Рс=0.2 pu Ос=0.1154 pu 
К,с=-2.875 К›с=0.2 
Кзс=0.2 Kyc=-2.681 
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3. MODELING OF INTERCONNECTED MICROGRIDS WITH LVAC NETWORK 

At steady state, the power transferred through the interconnection remains same at all conditions. 
Any change in demand or generation is taken care locally by the respective microgrid system. The system 
line diagram is shown in Figure 1. 


Рын ә. V12461, Fl 


PMSG Cos 


Figure 1. Line diagram of MG-1 and MG-2 connected through AC line 


The power balance equations are given as (1) (2) (3) and (4). 


PBG1+PWES=PL1+Ps (1) 
QBG1+QWES=QL1+QS (2) 
PBG2+PPVS+PG+Pr =PL2 (3) 
QBG2+QPVS+QG+Qr=PL2 (4) 


For small change in power (1) (2) (3) and (4) can be rewritten as (5) (6) (7) and (8). 


APBG1+APWES =APL1+APs (5) 
AQBG1+AQWES=AQLI1+AQs (6) 
APBG2+APPVS+APr=A PL2 (7) 
AQBG2+AQPVS+AQG-AQr =A PL2 (8) 


The change in active power generated and absorbed by load will result in the system frequency deviance 
which can be expressed in s domain as (9) and (10). 


AF 1(s)= = -[APBGI (s)+APWES (s)-APs(s)-APL1(s)] (9) 
AF2(s)= = =[APBG2(s)+APPVS (s)+APG(s)+APr(s)-APL2(s)] (10) 


Similarly, for reactive power the equation can be written as (11) and (12). 
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AV 1(s)=———[AQBG 1(s)+AQWES(s)—AQS(s)-AQL 1(s)] (11) 


1+sTv1 


AV2(s)=— [AQBG2(s)+AQPVS(s)+AQG(s)+AQr(s)-AQL2(s)] (12) 


1+sTv2 

Рвс, Pwes, Рруѕ, Pc, Ps, Р: and PL are active power of biogas, wind, solar, grid, line sending end 
power, line receiving end power and load. Овс, Qwzs, Оруз, Ос, Qs, О; and Qi are reactive power biogas, 
wind, solar, grid, line sending end power, line receiving end power and load. AF1 and AF2 are change in 
frequency of MG-1 and MG-2 respectively. Kp, Tp, Kv and T, are power system gain and time constant for 
change in active power and reactive power. 


3.1. Modeling of WES and PVS connected inverter 

The WES includes turbine together with PMSG which generates power having frequency and 
voltage variable in nature depending upon wind speed. The generator side AC-DC power electronics 
converter transmits the power of the PMSG to the load side DC-AC converter through a coupling capacitor in 
between the converters. The PVS comprises Photovoltaic Cell coupled with Buck Boost Converter. The PV 
side DC-DC converter transmits the power of the PV cell to the load side DC-AC converter through a 
coupling capacitor in between. The coupling capacitor helps in maintaining constant voltage across the 
convertors. Small signal model of the WES and PVS is established using power output equations from the 
Inverter model given as (13) and (14). 


piny = TUTTO) (13) 
QINV = im Cos (bv (14) 
For small variation in the output power of the Inverter (13) and (14) can be written as (15) and (16) 
APiv(s)=KinviAd(s)+Kinv2A0(S)+K inv3A V in(S)+K inva V(s) (15) 
AQwv(s)=K invsAdin(S)+K inveA@(S)+K inv7A Vin(S)+K invsA V(s) (16) 


The parameters/constants Pin, Vis.Zó and Xrwv are input power, inverter voltage and reactance of 
coupling transformer. Та is system time delay. In case of WES Tg=2.5s and for PVS Та=0.005$. Kinvi, Kin? -- 
-- Kinvg are constant of the above-mentioned equations, associated with state variables. The block diagram of 
the inverter with real and reactive power control is shown in Figure 2. The function of the proportional 
integral controllers is to establish Дб and AViny of the inverter such that APinv=AP’iny under steady state 
conditions. Data pertaining to WES and PVS are given in Table 6 and Table 8. 


AP "e(s) 
APin(s m — —» Kinv1=Kinv2 
Е 1 АРіпу(ѕ) 
FOPI3 я Kinv3 — Е 
l+sTe 
A6(s) 
AAS UA > Kinv4 
+ 
DCP) КаК 
AQinv(s) 
> Kinv8 > 
AV(s) >| FOPI4 > Kinv7 


Figure 2. Transfer function model of inverter 
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3.2. Modeling of biogas generator set 

The block diagram of biogas system having electronic speed governor is shown in Figure 3. Tac: to 
Твст are time constant of biogas engine mechanism. Data pertaining to ВСІ and BG2 are given in Table 5 
and Table 7. 

The regulatory valve supplies biogas to the engine. Valve position is regulated by the actuator which 
depends on signal from the speed governor. The signal is generated in response to frequency and voltage 
angle change. The change in active power output of the biogas-genset is given by (17) and (18). 


APBG (s) = ——AXVB(s) (17) 
Е (1+5ТВСА) 
AXVB(s) = s(1+sTBG5)(1+STBG6) AXAB(s) (18) 


The regulation in the actuator input signal is governed by the speed governor design and it is given by (19). 


(1+sTBG3) 
s(1+STBG1)(s2TBG1TBG2) 


AXAB(s) = APABG(s) (19) 


An IEEE type-I excitation control as shown in Figure 4 is considered neglecting saturation. The small 
deviance in voltage behind transient reactance, AE qB , is given by (20). 


(1 + sTgg¢)AE’ авс = [Kig ДЕғав (5) + KogcKAV(s)] (20) 
For small deviance in system voltage, the reactive power added Бу the biogas-genset is given by (21). 
ДОвс (S) = [КзвсАЕ'ав(5) + КавсДУ(5)] (21) 


Assuming that SG is operating at constant power factor, the reactive power change due to change in real 
power is given by (22). 


AQ'sc(s)=Kip APsa(s) (22) 


AP sc 


* 


1+5Твоз АР 1+5Твоз 4 1 
(1+5Тво + s/Tsc1T562) 5(1+5Тве (1+ 5Твоѕ (1+ Тве7) 


Electronic Soeed Governor Actuator Enaine/Generator 


Figure 3. Transfer function model of biogas-genset for real power control 


V(s) 
P 28 
x +)— Vret(s) Saturation Function 
<> r 
AV(s) "mu 
«і Кав <-> 1 
\ \ Y TS 
C» 1 + 5Тлв < ) Keep + ST ep 
V 4 d АЕгав (s) 
Amplifier Exciter 
Stabilizer 
SKrp 
1+sT rg 


Figure 4. Model of IEEE Type-I excitation control system for biogas-genset 
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3.3. Modeling of main grid 

In a grid connected system, the transformer with the help of online and offline tap changer is used to 
regulate the secondary winding voltage. Primary voltage of the transformer is considered to be constant. The 
slight alteration in active and reactive power of the grid can be written in state variables as (23) and (24). 


APo(s) = KicA@(s) + K2cAV(s) (23) 
АОбс(5) = K3cA@(s) + K4cAV(s) (24) 
The change in system bus voltage angle with change in the system frequency is given by (25). 


Ad EAE Q5) 
M 


3.4. AC network 

Figure 5 shows the model of AC network. The AC interconnection is a low voltage short ariel 
bunched cable to interconnect the two clusters that we are considering, 1.е., MG-1 and MG-2. Transfer 
function model AC Network in terms of power flow is given in Figure 6. Ps and Qs are sending end active 
and reactive power and Pr and Qr are receiving end active and reactive power. Being a short line, line 
efficiency is considered 9596 and power at receiving end is 100 kW, therefore, the AC Network have been 
designed considering line loss of 5 kW. The impedance of the line is given by ZZV-0.6232725.55*. 
Incremental change in power flow due to change in load or change in generation is given by (26), (27), (28), 
(29), (30), (31), (32) and (33). 


= vg Cos w—v1v2Cos(V +642) 


Ps 7 (26) 
Qs = Ир раза) (27) 
рү = eitaCosC ra vico (28) 
Py арза (29) 
APs = K1acAV1 + K2acAV2 + K3acAó12 (30) 
AQs = K4acAV1 + K5acAV2 + K6acAó12 (31) 
A Pr = K7acAV1 + K8acAV2 + K9acAó12 (32) 
AQr = K10acAV1 + K11acAV2 + K12acA612 (33) 
Where Klac, K2ac, ---- K12ac are constant of the above-mentioned equations, associated with state 


variables. The values of the constants are given in Table 2. 


3.5. Fractional order proportional integral controller 

In this paper, the simulation has been carried out using fractional order proportional integral 
controller. For FOPI FOMCON toolbox in MATLAB/Simulink is used. The improvement/tuning of control 
system performance in the time domain is carried out by minimizing error signal e(t) using integral square 
error (ISE). The fractional order PI controller transfer function is given by (34). 


G(s)a = Kp +=; (34) 


Where a is non-integer. In case of FOPI three variables i.e. Kp, Ki and a need to be tuned instead of two 
variables in case of conventional PI Controller. The values of Kp, Ki and а are given at Table 1. 
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avi 


0.754 Mhetat2 K7ac 


Figure 5. Model of AC network 


4. SYSTEM SIMULATION AND RESULTS 

The transfer function model of the microgrid system MG-1 and MG-2 using PVS-biogas and PVS- 
biogas-limited grid is shown in Figures 7(a) and (b) respectively. In this paper only one typical example is 
considered i.e. in both the Clusters MG-1 and MG-2, step load is increased by 1%. Further, it is assumed that 
input power to WES in MG-1 increases by 1% and input solar power to PVS decreases by 1%. Figure 7(c) 
shows dynamic response of the system (MG-1) due to 1% increase in active and reactive power of load and 
1% increase in wind input power, time (second) versus change in voltage. 

In MG-1 due to step increase in load and increase in input wind power by 1%, the wind active and 
reactive power output increases gradually along with decrease in output of biogas. It takes around 9s for 
output power to reach steady state in MG-1 in Figures 7(d) and (e). This is due to the slow response of wind 
turbine. However, frequency and voltage angle reach steady state within 5s in Figures 7(a) and (b). The 
voltage reaches steady state within 60msec in Figure 7(c). Figures 8(a), (b), (c) shows dynamic response of 
the system (MG-2) due to 1% increase in active and reactive power of load and 1% decrease in solar input 
power, time (second) versus change in frequency, change in voltage angle and change in voltage. 

In MG-2 due to step increase in load by 196 and decrease in solar power, active power output of 
PVS inverter decreases immediately in Figure 8(d). To compensate, initially power is drawn from the grid 
which slowly decreases as the output power from the biogas increases. The steady state of active power is 
achieved within 5 sec. Also, the active power drawing from the grid becomes zero. However, the case of 
reactive power in MG-2 is different in Figure 8(e). Even though the solar input power decreases the coupling 
capacitor attached with PVS inverter provides reactive power initially and finally taken up by biogas 
reaching steady state within 5s. As seen from the graph the frequency, Figures 7(a) and 8(a) and voltage 
angle, Figures 7(b) and 8(c) settles within 5 sec and voltage, Figures 7(c) and 8(c) settles within 60msec in 
both the microgrid as anticipated. Further, from the frequency graph of both the microgrid it can be 
interpreted that frequency oscillation is more in MG-2 as compared to MG-1; this is because of presence 
WES inertia in MG-1 which supports the system working as shock absorber. Oscillation in active and 
reactive power PVS and Biogas in MG-2 can be interpreted as system interaction with the grid. Coupling 
capacitor attached with the inverter provide the reactive power initially through the stored energy therefore, 
voltage profile is maintained quickly. As desired no extra power is drawn from the grid as well as AC 
network. From the above, it is opined that controllers are tuned to overcome the small disturbance in the 
system. 


Microgrids dynamic stability interconnected through low voltage AC network (Vinit Kumar Singh) 


AF1 (pu) 


334 o ISSN: 2252-8792 


Microgrid MG-1 Microgrid MG-2 
Ad 


A 


52 
: | APaw(s) Va APaw(s 
А Electronic Bio Gas i vA Electronic Bio Gas Fai 
Speed — | > Actuator & Speed [> Actuator & 
Govemor Engine Kyevs Governor Engine 


РИ РІЗ PIB PIT APszox(s) 


AR АР рг(5) 
m Keo PVS № 
S | AF2 1+8, Aon 
ДР: £ 


GRID 
А 


V2 


Excitation 
Control of SG 


Excitation | 
Control of SG 


pem 


Figure 6. Transfer function model of: (a) MG-1, (b) MG-2 
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Figure 7. Dynamic response of the system (MG-1) due to 1% increase in active and reactive power of load 

and 1% increase in wind input power, time (second) versus: (a) Change in frequency, (b) Change in voltage 

angle, (c) Change in voltage, (d) Change in active power of ВСІ and WES, (e) Change in reactive power of 
ВОТ and WES 


Int J Appl Power Eng, Vol. 10, No. 4, December 2021: 326 — 336 


Int J Appl Power Eng ISSN: 2252-8792 o 335 


a 1 
2 0 05 
5 * 
e 9 - E 8 о 
N ES a 
4-2 © 2 
4 -0.5 
-4 < 
4 
-6 
-8 А 
0 1 2 3 4 5 6 ò Я 5 3 4 5 6 о 0.01 0.02 0.03 0.04 0.05 0.06 007 
Time (seconds) Time (secs) Time (seconds) 
(a) (b) (с) 
103 
8r. 
n 
era 
34 fà APpv 3 
gyi re 2 
B ai ===> 3Р9 Я 
uy 
о | NEP SÉ 
2 
0 1 2 3 4 5 
Time (seconds| 
( ) Time (seconds) 
(d) (e) 


Figure 8. Dynamic response of the system (MG-2) due to 146 increase in active and reactive power of load 
and 1% decrease in solar input power, time (second) versus: (a) Change in frequency, (b) Change in voltage 
angle, (c) Change in voltage, (d) Change in active power of BG2, PVS and grid, (e) Change in reactive power 
of BG2, PVS and grid 


5. CONCLUSION 

Interconnected microgrid have modeled using wind energy system and photovoltaic system 
accompanied with biogas based system. One of the microgrid is also connected to main grid. The complete 
system is tested with 1% step change in input power to renewable energy sources and load. As per the result 
the system becomes stable i.e. frequency and angle settles within 5s and voltage within 60 ms. As desired no 
extra power is drawn from the main grid as well ac interconnecting line. Each microgrid gets support from its 
own system to override the disturbance caused by change in load and input power to renewable sources. 
From the study it can be interpretated that if the individual control of the cluster is achieved then 
interconnection of microgrids could be solution for valid optimization of resources. The benefit of 
ineterconnection is that at time contigency the power transfer can get reversed as per the requirement. Also 
the line supports initially when there is step change in load or input to the generating resouces. 
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